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SUMMARY 


This is an investigation of the equations of motion and 
physical parameters involved in stabilizing the initial flight 
of a vertically launched rocket by means of a booster rocket 
pin-connected below the main rocket. The system is designed 
to stabilize the flight in its early stage before the aerodynamic 
control surfaces become effective. Stability of the system is 
dependent on the pendulum action of the booster rocket. 

The equations of midi ee derived from Lagrange's 
generalized momentum equation. The differential equations 
thus obtained were not solved but were tested for stability by 
means of Routh's stability criteria. The ratio of the mass of 
the main rocket, Mj, to the mass of the booster rocket, Мә, 

My М; 
was investigated for the two values 77~ = 1.5 andy = 71:98. 

The system involving a mass сано ту, = 7.75 was found 
to be unstable under all conditions. However, the system in- 
volving a mass ratio * = 1,5 was determined to be stable in 
the range 1.62 44 < n x 1010, where g is defined as the ratio 
of the distance 4 from the center of gravity of the booster M, 
to the pin connecting the strut to the main rocket Мү, divided 


by the radius of gyration, Ку, of the booster м,. In this range, 


for any given value of Æ , stability was uniquely determined by 

















ШЕ» э o i 1 
hough the system was found to be theoretically stable 


for the mass ratio 7 = 1.5, the ratio 2 turned out to be of 
42 


such great magnitude as to make the system entirely impracti- 


calfor this particular mass ratio, 
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I, INTRODUCTION 


When a rocket or guided missile is initially launched there 
is a short period when its flight is unstable due to the fact that at 
low velocities the aerodynamic control surfaces are ineffective. 

Various procedures have been adopted to stabilize this 
initial part of the flight. The most obvious and most successful 
of these has been the use of guide rails, or launching towers which 
hold the missile on its course until it has attained a velocity at 
which the control surfaces become effective. But this method 
requires a bulky launcher, and special care is necessary to main- 
tain the rails in proper alignment. 

A modification of the guide rail system is the ''short- 
length" launcher which is really a guide rail whose length has 
been reduced to a minimum by greatly increasing the initial ac- 
celeration and thus the velocity of the missile. This increased 
acceleration presents a great many problems in component design 
to resist the tremendous accelerative forces. 

During World War II the Germans were fairly successful 
with a unique approach to the problem in their well-known "V-Z2" 
rocket. (Reference 1), Тһе "У-2" was launched from a near- 
vertical position with no external restrictions. Instead it em- 
ployed four carbon vanes mounted in the jet stream and activated 


by the auto-pilot to maintain the missile on its proper flight path. 








32. 
The great difficulty here was that the vanes steadily burned out 
with a consequent diminishing of control. 

The subject of this thesis is a third means of initial flight 
stabilization based on the concept of the action of a pendulum and 
independent of both external guides and internal control mechan- 
isms. 

If the missile booster is constructed as a separate unit 
and attached to the main rocket by a pin connection a certain uis- 
tance below its own center of gravity it can be seen from Figure 
l that any motion of the main rocket will impart a proportionate 
motion to the booster. Specifically, if the main rocket rotates 
about its center of gravity in a clockwise direction, the connec- 
ting strut will cause the booster to rotate also about its own center 
of gravity in a clockwise direction. Thus, if the main rocket 
turns off course in a clockwise direction the booster immediately 
turns in a clockwise direction also, that is, the booster tends to 
re-align its thrust with the axis of the main rocket. If, bya 
proper choice of moment arms, the booster can be made to over- 
shoot this position of thrust alignment slightly the thrust will have 
a horizontal component which acts on the strut to the main rocket 
and causes the main rocket to turn about its center of gravity in 
& counter-clockwise direction. Thus, the booster responds to an 


error and its response tends to counteract the error. Sucha sys- 








aĝa 
tem should produce a sinusoidal flight path if the proper moment 
arms and damping forces were applied. 

In this analysis frictional and aerodynamic forces were ne- 
glected to simplify the equations. The only damping force consid- 
ered (other than the inertia forces of the system) was the jet damp- 
ing force or the resistance of a jet stream to rotation. The other 
external forces considered were the weights of the main rocket 
and the booster, and the thrust of the booster. (The motor of the 
main rocket is considered not to be operating in this analysis). 
Taking into account only these forces the writer derived the equa- 
tions of motion of the system in two dimensions by means of La- 


grange's generalized momentum equation. Because of the com- 


plicated nature of these equations no solution of them was attempted. 


Instead the conditions for stability were investigated by means of 
Routh's stability criteria, 

A hypothetical system consisting of a five-second booster 
pin-connected below a ''V-2" rocket was first investigated. For 
this system the ratio of the mass of the main rocket My, to the 
mass of the booster M,, is = = 7.75. To reduce the work of 
computation the lateral motion of the center of gravity of the sys- 
tem was set equal to zero and only the vertical motion of the 


center of gravity and the rotation of the main rocket and the 


booster were considered. 












| = Ў 
ity was dete зай it was of such a sensitive nature 





is to make impractical further inve stigation for this mass ratio 


| the general case with lateral motion taken into account. 











Pocket, the weight, Mog, of the booster, the thrust, F, of the 


booster, and the jet damping force, D, of the booster. These 


forces are indicated in Figure 1. 


Figure l. External Forces on Main Rocket and Booster. 
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The jet damping force, D, is a force which resists the ro- 
tation of the booster about its own center of gravity. If the booster 


is of mass Mo, then the rate 


77 


of flow through the nozzle may 


ам, З 
be written as a ‚or м». As 


the booster rotates about its cen- 2 
ter of gravity the tangential ve- 


locity at the nozzle exit will be | 
f 
"E or 1, 4, . Therefore GL 


the damping force, D, will be 


equal to the rate of change of momentum, or м; 20:3 q, ). This 





M 
expression has the dimensions т 1-2 ‚or Эр . 


Frictional forces and aerodynamic forces are neglected in 
this analysis. 

Both the main rocket, Mi, and the booster, м. are ав- 
signed three degrees of freedom, namely, lateral, vertical, and 
rotational motion. Note, also, that only the two dimensional case 
is being treated. The coordinates of masses My and M, are indi- 


2 


cated in Figure 2. The connecting strut between M, and M, is as- 


1 
sumed to be of zero mass and of infinite stiffness for purposes of 
this analysis, and consequently the coordinates of the center of 


gravity of the system are indicated on the straight line joining 


Mı and M5 in Figure 2. 











Figure 2, Schematic Diagram of Coordinate System and 
Dimensions. 


эӊ 
It is assumed that the propulsion system of the main rocket, 
Mj, will not be operating during the initial stages of launching and 
consequently м, is a constant for this analysis. The booster, м. 
will of course be burning with a consequent change of mass during 
launching, and the rate of change of mass M5, say M5, will be 


considered, as will be the change of position relative to M, and M5 


1 
of the center of gravity of the system. 

The first step of the analysis is to eliminate the extraneous 
coordinates 45, 46, 47» and 98 in order to write the equations of 
motion in terms of the rectilinear coordinates of the center of grav- 
ity of the system, чү, and qz, and the angular coordinates q} and 44 
of M, and Bü, about vertical axes through their respective centers 
of gravity. (Since only the onset of instability is of interest, the 
analysis will be —— to consideration of small variations of 
the angles q4 and q4.) Once this step has been completed, La- 


grange's generalized momentum equation will be used to deter- 


mine the equations of motion. 


ELIMINATION OF EXTRANZOUS COORDINATES 
The coordinates q5, чє, ау, апа Gg are always redundant 
because of the rigid links in the system and the known position of 
the center of gravity of the two masses. Therefore these extran- 


eous coordinates may be eliminated by considering the constraints 








elaborated below. 

The first relation between the coordinates of the individual 
centers of gravity and the center of gravity of the system may be 
found by taking moments about the center of gravity of the system 


indicated in Figure 3. 





PIN 


Figure 3. Relationship of My and M, to the center of gravity of 
the system. 


Thus 77 = Ag RI =0 


or ИФ) = МО) (1) 








ES 
Furthermore the horizontal and vertical separations of 
the centers of gravity of M, and M, in Figure 4 are fixed by the 
length of the line, 7, joining the two, and by the angle ,/q, - 2 


which 2 makes with the vertical. 











Figure 4. Relationship of M; and м. to the connecting strut. 


In Figure 4 the vertical separation is 


je fa = ЄС ees (p-d) 


or 4 = AED (2) 


and the horizontal separation is 
hr 75T 


EU 22% т 28054) (3) 


Me 
A fourth equation follows from a consideration of the 


similar triangles indicated in Figure 3. 


Е 5 с/т 

fe “fi AU 

— г = 24274 0 С3н 
Mge- g? ^h (Fe) 


and, since the denominators of both sides of this equation are 


equal from Equation (1), it may be written 


ARA) = Me) 5 

Equations (1), (2), (3), and (4) provide four linear rela- 
tions between the coordinates 4,.4,.4,. and q,, andq, ,q,, 
9, , and q, >» Consequently the former may be expressed in terms 
of the latter and thereby eliminated from' subsequent equations. 


Solving these equations the following expressions are obtained: 


From Equation (2) yo em L, «03 (7,-4/ + fa (5) 
From Equation (3) ғ 4y +2 Ga -4) + 7” 
ог | Boc ge rm 07 4) (6) 


Substituting Equation (6) into Equation (1) 








cles 


AE) 3 ^R 6 »,(-4)2 
leo зо ой Ser * 144 2 0-7) 





f* = y + Ps A vg, d) (7) 


Pt IM 
Substituting Zquation (5) into Zquation (4) 
ED) X cos (7-4) + fe - ff = TA 


Beta) = MR +17 -^ US Ly cos (5-5) 


ж 


ft Е л 3 ^ot ^ 4 aa G4) (8) 





Substituting Equation (7) into Equation (6) 


ff = 7 + pe “С э» (ў,-4) 39 aim (Ga - 4) 





Rot p e 1) 


2” — E (9) 
From Equation (2) 
Ж = 22 es» (fs - 7) Tr Ж 


Substituting Equation (8) into this equation 
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ГЛЕ ЕИ 
5-5 йыры Cy 
Ғы uc 77 (10) 


These expressions for the extraneous coordinates in Equa- 
tions (7), (8), (9), and (10) involve, in addition to the principal co- 
ordinates, the variables -Z,, the distance between the centers of 
gravity of the two components, and Л , the angle between the strut 
4, , and the line 4, joining the centers of gravity of the two 
components. These two variables will now also be expressed in 
terms of the principal coordinates. 


In Figure 4 by the law of cosines 


-% ш ВЕС 2121-2214 eos Gr J=? (11) 


and by the law of sines 


ыг . 224 


MU - Ж; 
epee ©, 
Jag = LA sia Fe Je) 


D a f e 293A 


whence fc = f = A NA (12) 
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These relations are unfortunately complicated. There- 
fore, to simplify them for later use, the fact that the variations 
inq, andq, are small will be taken into account. Assuming 


these angles sufficiently small such that terms of order q 


may be neglected, the following approximations can be made 


310 * E and equus “ / 


From Equation (11) 


Gb =\ 474-214 = L La 
From Equation (12) 
(4-4) = ғ 
6-4) 
0) = д4 - Ф т 


and since ty 22 К e 





fe rp) 


o 





= = — - La = Z - 
Соз (% = г) 
Thus the appropriate approximations are 


co / fa -4/ 2217 


ES Ж, 4 
sin (3, -d) 22 2-4 fr > LA f“ 


— А 


9 


(13) 
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Substitution of the results of Equation (13) into Equations 
(7). (8), (9), and (10) yields the following approximate relations 


between the extraneous and the principal coordinates 


e 


+ 





А, 
227 мем СОЕ 00) (4) 


ДЕ 


RR (A) S 





С Oe eu (16) 
2.52 з Ба (17) 


To derive the equations of motion by the method of La- 
grange the expression for the kinetic energy, T, must first be 


found for use in the Lagrangian equation 


IT OT 
Fa J (18) 


where E^ is the "generalized force" and 4, is ЕЖ (Reference 2). 
The generalized force 2. is employed rather than the potential 
energy because this is à non-conservative system. 

If I; denotes the moment of inertia of the mass M, about 
its own center of gravity, the kinetic energy of translation and 


rotation for the system may be written 
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T= tls ther ER) = af + Je) (19) 


DERIVATION OF EQUATION OF MOTION 


FOR COORDINATE а, 


In applying Equation (18) to the direction of the coordinate 


qj» the value of -2£ is first computed using the kinetic energy, 


ГА 


T, of Equation (19). 
27 =~ Hg Be - m. 
% Ё 2f, f of, 

The extraneous coordinates q, and 4, шау be eliminated from 


this expression by differentiating Equations (14) and (16) respec- 


tively. From Equation (14) 








| 046. ж 1 Е м н 
57 pom (Spe Afi) Bey, A p np) 
ЖЕ, 
2 

and from Equation (16) 


К ДОЛ ЕЕ хэ, ж. 








бт 
2А 2и 
of 
Therefore 
Н) а уч 
— AE o fee ORAR) y AN 





Al t IA n uu p P eg 


(4-74 








2132 


ташы — ыг, VR) ^ Gi —— 5927 


The first term in the Lagrangian expression is then 





ayer] Ша ЛЯ). адар ан) 


Af! AF (те) 


ҰЛТ» / 2 de (Af 4f] - 2255, (49 -4%/ 


(447% )* 





ML xc yy Agr) 





мМ 


22 + 202 Mth (bp n 
or er) = (raji Af RÀ) 





er (Ly -A g) + E E (20) 


By inspection of Equation (19), the second term of the La- 
grangian expression is 


27 = (21) 
27 


The generalized force 22, is most easily computed by con- 
sidering the work done by the external forces as the coordinate q, 
is varied, all other coordinates remaining fixed. 

Work; 
Afe 


In particular, the force & becomes, referring to Figure 2, 


25 2 Work - [E sandy + D cos Уу) A 2 
p ^ f. ag 


In general, Work, - X "E . OR ec 


” 
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But within the present approximations 


227 2 and 203 Ууж / 


hence &, = 22 + 2 (22) 


Substituting the results of Equations (20), (21), апа (22) 


into Equation (18), the equation of motion in the q, direction be- 


comes 





е) + hile Np Af) - fg Ap) 


+ e = С + 2 (23) 


DERIVATION OF EQUATION OF MOTION 


FOR COORDINATE q2 


The equation of motion in the 4, direction is obtained in 


an identical manner. 


From Equation (19) 
IT Mia „ ^45 M 
and, differentiating Equation (15) 
2 Е Ж mn (4- 4) 
Consequently Ae = д 
of 


Similarly, differentiating Equation (17) 
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2 = Ф + fe. (A) 





so that 
fe = / 
л 
Therefore 
2 : Mt pts 
28, = (nern) [fa + (7% 2227 (4M 


and the first term of the Lagrangian expression becomes 


EM 2 MA ero А [LAI] A др 4-47 


or 


ap) = areca та 


Furthermore = O (25) 





27” 
fa 
The generalized force 2, is, referring to Figure 2, 


Ж -Жоғк . 420%-74/9 «б<ез%; - Din paul аук 
— Б 5 


and since 


3/8 fa E fr and сэг fy = 


A = -Co атг? (26) 








cA 
Substituting the results of £quations (24), (25), and (26) 
into Equation (18), the equation of motion in the а, direction be- 


comes 


(ten) ru bho ag =F - gern) - Фу, (27) 


DERIVATION OF EQUATION CF MOTION 


FOR COORDINATE ч, 


The first term of the Lagrangian expression is obtained in 
an identical manner to that employed for coordinates q, andq, . 


From Equation (19) | 

— = L, fr + 227 of: 7 ^t f, 22- 

{з 2% "e 

and, differentiating Equation (14) 

АҒАР АУА АР а ps tk) 
so that 

Ф - 4% 
Similarly, differentiating Equation (16) 


лт NS 


JG OR GD + Ge 


so that 








3213 


Therefore 


% - 2440: 7 127277 HERG „ЖИ 260) 


ж /2- 5% Np + pty. Up p ae) 


or 





a = Lp + ¿LE L(Lfo Af) 


and the first term of the Lagrangian expression becomes 


AGI “Lf, E ФЕ) * A Ug AR) (28) 


From Equation (19) 


and, differentiating Equation (14) 


fo = fit at Ud Rr Ug 
so that (30) 


ef? B (7 +12) 


Similarly, differentiating Equation (16) 


ARALAR тєр p-ga) 


so that (31) 


ae = ^t P, 4 








But now the values of fe and 2: cannot be evaluated to 
the first order from the linearized expressions for 4, and q, 
in <quations (15) and (17) respectively because the q, term has 


vanished in the process of linearizing these equations, Therefore 


the differentiation must be carried out upon the exact expressions, 


Equations (10) and (8), and the results then linearized. 


From Equation (10) 


ис зу “ж = cos O where © 22 


and the time derivative is 


(32) 








==: : 5 ж E 32223... о 
je Ё а, бөз. А беө. А, My cas 


Therefore 


2h 2 1% LA omod б дене 5. _a_ (mo LS. 








2% Ж», ZA Ме 2f; м Ма 77 
+n 24 соз © Al ад dcos O ^a A E oO E 
rte Of” ар "ойлар Зу 

^7, ^7 соз 
2,6 (33) 


4 —— 
АЛАДЫ ZE 


The value of 2% follows from Equation (11) 
^f 


4, = 7 A = г 2 Е eos (fy - ә) 








8573 


and upon differentiation 


2” = 5 228 sia (H-B) 
2% Ly 


which, by retaining only first order terms may be written 


P w ALTE 


Ya Ly 
Similarly (34) 


24 _ 2 _ Atlante pal 
2” 27 


The various derivatives of © must be computed from 





Equation (12) 


o = pd = pi- oa LA sarph 


Thus 


й / боз 2) 


EU Te za 


do 26 Ly cospe) -Lh за eds) 
4 qu Ay 


and approximately 





— 2 210120 12 (35) 
273 Ту ЖБ “ S 








Now 


25/29 = Соз oO au 
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and substituting the results of Equation (35) and linearizing 


Dane „Lesen, ш (36) 
29, p Ay 


Similarly, from the results of Zquations (35) and (13) 


2038 317029 ~ ФА £L 22 X nA p) (37) 
ж” Ay La? 


^ 75 


Substituting Equations (34), (35), (36), and (37) into Zqua- 


tion (33) and retaining only first order terms 


Р ЖЕ & -26%--ӨЖЯ яа Lb (tui) 
ME. = _ u AT V) — — 


2f; ж ty жж” 


ж”ь ELA - Ф) = ALIA L Fa -6 
C (n4) p (TAIAN УЛ 


or 


- 


— ж LA hf) мм E Ру (38) 


29, ^ 5 (0) 27 е, 


To find > the same procedure is followed. Thus, from 
? 


Equation (8) 





=e — mm = gz- 
e Ge NO. Ж = where © f 


and the time derivative is 


n 25 Son d EPA E. pp mr cos © (39) 








Therefore 
le ^h 248 ¿na E y ty ¿IO Ф 
A ^T 25 


IB = 
2f /Жжеғ?а “ 


А, EL. Ipe IS 32 7, a4 соз © 
IT vw LT fs 

















ar м + ТА 3 
17, 9^ PITE AL IT 2% cos 9 yA BA eos O (40) 
CAU п От)? дт (ra)? af 


Substituting Equations (34), (35), (36), and (37) into Equation (39) 


and retaining only first order terms 


Be 2 Lip tip) + qe Se (fr fr) 








22 М оа 
20 E = — Lips haze) 


27 -- £ 2 ^ Mr er ^ — 7 2-4 
263 wx 6 ав Е 





The second term of the Lagrangian can now be written by 


substituting Equations (30), (31), (32), (38), (39) and (41) into cqua- 
tion (30) 


т 
T СТ 17) 


( MR ea га A ify th Jr) 4 28 en 203 “би ан 242%, 47 


+ fe ae OP Ad fe Op ӨМ era A] + 


moa 








— 





AB ERA th cta 





+ Elfe a С 


or 





PC X 
ж” — С гау AE fo 2 
— С (42) 


The generalized force 2, is 
Z, = Work, 


afo 


Therefore, referring to Figure 5, 


AK IA 2 77 
A E иас Ж 2 Al Arh) FM 


ел Lj- latt) epg] 


+ Don "m (Arh) sin gaf 


2222 ЕС a) 








DI 
Substituting Equations (8), (9), and (10) to eliminate the extraneous 


coordinates Ч, » 8,» and 4,» and collecting terms 
52, ССРИ zem zen 4 3-€09- (hh) son pu] 


+(Fees 2 - 22,9) Zo 5 L cas © -2-0) ees fy] 





pe mew 


Figure 5. Schematic diagram of external forces acting on system. 


Carrying out the indicated differentiation with respect бо с, 


22, = (Fain Pos zu) — o (-6 “2” == N 


«(A msg, - 2 310 Pr) [- Сл 6 — соз © 7] 








эс 
=£0- 


Substituting Squations (34), (36), and (37) into this expression, 
and using the approximations sn f. — and vore E o! a 


becomes 





= @ +o) f- utu * AA. 27 = A pel n 


(Ela AR ше || 


Finally, retaining only first order terms, the expression for the 


generalized force <2, reduces to 


— AAA (43) 


M, AAN 
Therefore, by substituting the results of Equations (28), 


(42), and (43) into Equation (18), the equation of motion in the q, 


direction becomes 


О оО 00-40. яны <%-а)) 








— x x жағ, Ip) e LP 45 ў 
7-2 404 (0-4) = AT 6 Ге) +27 (44) 
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DERIVATION CF EQUATION OF MOTION 


FOR COORDINATE q4 


The first term of the Lagrangian expression is obtained in 
an identical manner to that employed for coordinates q, .9,+ 


and q,. 


From Equation (19) 


al = af + CREE» Бог. 
2f, f 224 mu 


and, differentiating Equation (14) 


f 





Pe ==? AAA) AR Ap Ap) 


so that 


27 -. жб 
Ж 14 +772 


Similarly, differentiating Equation (16) 


Be ж Й 60) + Lp Ap) 


so that 


Therefore, 





r . : nts 
E = La fa ж ELR A eA (“%- 57:22 Ааа а ЈЕ ] 





— men -4 f) a ya ES APN AZ] 


or 


2роэ up - Ф 0-4) 
Y 
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and the first term of the Lagrangian expression becomes 


SE EL bs RA RAR 


IN 
% 


-a 4 fe) (45) 


From Equation (19) 


Zola) pgp ag) © 


and, differentiating Equation (14) 





u ler.) 
so that (47) 


2 - _ ола 2 
ж 7 AP 


Similarly, differentiating Zquation (16) 





м 
# = d 221200 -A fa) d mar R MR) 
so that (48) 
b=. LEN 
E c pue 
24 әр : 
Here again the values of ЭР and 2% cannot be evalu 
ated to the first order from the linearized expressions for q апа 


q, in Equation (15) and (17) respectively because the q, term has 


vanished in the process of linearizing these equations. Therefore 


= 
the differentiation must be carried out upon the exact expressions, 
Equations (10) and (8), and the results then linearized. 


The time derivative of q- , defined by Equation (10), has 


already been found from Equation (32) 


2 = № - ясы As, (Ace 9)ó “хэ A Apes Oy Ж ОЮ) 4.932) 








Therefore 
A = e edt (arm o/ = a a 23140. Ò 
27% “627 2 Jy ^7 + 1, уе 


A e 78 Ma 44% сызсо 
МАТА am 2 * ^7 2fv 





Ma $ darb p Lath 24,499 gr (7 сез 
2 2f» M Pv SL p^ Qum)? дуу (49) 








2% : - 
The value of E follows from Equation (11) 


27 = A e EAS eos (fy fs) 


and upon differentiation 
2-4 = ACA 3/7 (Px -f3) 
279 ty 
which, by retaining only first order terms, may be written 


2242 x= BGA Tv - fs) | 


ae 4 








-32. 





Similarly 
oh E 2 20-44 (уугуз)7 СУУ - fa) (50) 
724 27 


The various derivatives of © must be computed from 


Equation (12) 


o pr o uv 


Thus 


PR pr lA 
Loe 


д Joa Fa 2) 7 > 
Y — a 22 Y 
^ 





r 
20 So 222 Arcos pega) -. s ) 
27” ANA, se Gr) A” 
and approximately 
24 12. — — (51) 
SE 24 2087 Ar 
Now 
2 зә ee = Los со eae 


“, of 


and, substituting the results of Equation (51) and linearizing, 
» 32-4 = I cos O ~ E (5 2) 


Similarly, from the results of Equations (51) and (13), 


Pese _ — = Ä 2; т La lata po) 53 
Av > 757 % 2Ф Жа (53) 
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Substituting Equations (50), (51), (52), and (53) into Equa- 


tion (49) and retaining only first order terms 








"utr Ze AR) + a ERE 


MIR LA Gazga) * Stn „4 (£ ga -A fr) 
> 


t Orern)" Jj? 27% 


or 


SR ALP (** ^7. 


To find = the same procedure is followed. The time 


Зу ote. A fe OA) e Vau ӨСТІ (54) 


derivative of q,, defined by Equation (8), has already been found 


from Equation (39): 








» а - 5 > 3 
jr = 48 it 27 44 (sina)ó - 22 are (39) 


ж»4?, 


Therefore 





JE. M IH (no) 6 “Ж СЫ уос ж? DE 
“Д7 лу — Жа O + 2 My om 2% 


“ - 2-4?» Же. 


— - 4 Iso 
APA yc А Ae 


AL AT 26 2 ^L ^. 9 eos S 
"@ mpm 92 ойы 9% z 
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Substituting Cquations (50), (51), (52), and (53) into Equation (55) 


and retaining only first order terms 





Pfa AT. ; А 09-Х 
ЯВ- BEA) = ль, Sole fo 
ET чш lu 


or 


202004) » E 4-4) (56 





огт м. 
The second term of the Lagrangian can now be written by 


substituting Equations (32), (39), (47), (54), and (56) into Equa- 











tion (46) 
dio PME UR A SA] 
А д ДЛ MÀ e EE Ups A] 
«Да, ы 5 + er Ad 
+ Agi а “бие 
- 
2-- а 47 - en A (App) 


— A RO A) (57) 


1%, 1%, > 
еи af + ЖЖ ГР 











sabe 
The generalized force A, is 
«6, = Work, 
4% 
Therefore, referring to Figure 5, 
Re = (Ben s ege) de [y (do) en jo] 
„(оез - D anpi) -2 5; Ih - 1222) 5922 
AI LE - жу Ж 
2” % БХ 2% 


Substituting Equations (8), (9), and (10) into this expression to 


eliminate the extraneous coordinates 49, .4,. апа 4, 


20 s Pas da o 4 a 2 Ay 200 8 - (he by) 2:0 pa] 





+ (43h, - Pf) LR - Ы 42-6 «з9-(4-.%) cos fol 





vA gf utut cu Ж МЕ ы 


and collecting terms and performing the differentiation 


= (©з; r Deor prf- ra Suter E — — ae 705202 





v (A955, - Ds. af) [- wA deum, cos 924)» (Lek) siap] 


Substituting Equations (13), (34), (36), and (37) into this expression, 


and using the approximations aap = and cos xl, I, becomes 
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-/ 25 (-4/ -(/4--4/ 


(7 tl ja (LU s Le Heep] 


Finally, retaining only first order terms, the expression for the 


generalized force 2) reduces to 


Ay = Dy M ou] (58) 


Therefore, by substituting the results of Equations (45), 
(57), and (58) into Cquation (18), the equation of motion in the q, 


direction becomes 


AE RA - AU HA) fe Ap Ap 





+ fe fe + AEB, 4-4) - FO RÀ 





- е, 45 64-5) - -2 44 4--57 (59) 


Thus Equations (23), (27), (44) and (59) are the four La- 
grangian equations of motion corresponding to the four principal 


coordinates of the system. 
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ANALYSIS OF THE EXTERNAL FORCES AND REDUCTION 


OF THE EQUATIONS TO DIMENSIONLESS FORM 


In carrying out the stability analysis the state of vertical 
motion corresponding to the q, direction will be considered con- 
stant during the perturbation. Thatis, the values of 4, and q, 
will be fixed, and their magnitudes will be determined by the 
time-instant at which the stability of the system is investigated. 
Therefore the equation of motion in the q, direction, Equation 
(27), will not be considered further in this analysis. 

The thrust force, F, of the booster may be expressed in 
terms of the mass rate of flow, M, ‚ and the specific impulse of 


the fuel, I, , as 


"e c - 1% q La 


The negative sign is necessary because M, is numerically a 
negative quantity as it occurs in the Lagrangian derivation. Thus 
the thrust, F, will become a positive quantity when evaluated nu- 
merically in the next section. 

Also, it was found in the first part of the derivation that 
the jet damping force, D, is proportional to the quantity EA Ja 
In order to make this a numerically positive quantity, it will be 


written as 


DP = mA 








The moment of inertia, I; , can be written as the product 


, And the square of the radius of gyration, k, , 
4 


of the mass, M, 


that is 
L = G4 


The time derivative of the moment of inertia, I., is 


Ж”; — . However, the second term is 





227 66), 
Yi ot 


small compared to the first term and will be neglected in this 


analysis. Therefore the following approximation can be made 


Making these four substitutions in the equations of motion 


in the q, , q, and q, directions the following equations are ob- 
tained from Equations (23), (44), and (59), respectively. 


From Zquation (23) 


Mali + rash “fe an 
(60) 


2222-1200 — +492.) я = о 


LUT + Ah 





From Zquation (44) 


EL] AAA AL le T со” UA T 





- = 
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А 272 Ж — айбаты E 47 <2) % 
AZ A fe ее. fit ЈА 


ж MA - 
+ (Ge: (> ae (61) 





From Equation (59) 








малма d ё“ IEAI, 2) SS 
— — ^7, — ⸗ E 
(242074 ПАИ Ея 4) fe 








ААА, Иа AA) fp = 2 vem 


To write these last three equations in dimensionless form 
introduce the dimensionless coordinates Q, , C, , and Q, and the 


characteristic units of time z lo velocity = I,,g, and acceleration 


ш Б» 
Thus 
“ = 5 Giz Ze 2 ш 
zs #9 5 
$-5 $5 Z =Z, p ee 


& 
u 
N 
4% 
ң 
N 
N 
5. 
& 
і 
N 
E 
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Making the above substitutions into Equation (60) and divid- 
ing the resulting equation by Mg, the following dimensionless equa- 


tion is obtained 


^l^ д ж Ma Lap ё, A 
bre Bld + Boa AA яг? л 











„Сла. Xx 26 O a AZ A. 4 
10, * ^f. => 4.4% >? Ме 7274 А 


(НЕТА a 
21220122 = 2 (64) 


Similarly, substituting the values of Equation (63) into 


Me Me 


Equation (61) and multiplying the resulting equation by 
Mgt 


the second dimensionless equation is obtained 


“gen Sl? ARA FAA 


wel 




















lar ke LLG ^7, 2те 475). Г, 3 58, 4 
а 22 285 (768) — 274 


мем 9 А А Ж, Л, “(Ж 


on 


+ Ar = 4 mo 








— 2 7 GE 22:76 2 - Le lir] a= (65) 


ж” 227 ^^, 2,99 4, Mote 
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Finally, substituting the values of Equation (63) into Equa- 


tion (62) and multiplying the resulting equation by — the 


third dimensionless equation is obtained 











т 22- д a А» 2 £L 2 
мж” 4g ^ 4 2 


226 thie k AAG pm [tr] А. .4 6 о, 
ж Я.А 27 64-4 7 


22222 2 (fr Zr 2 4/8 2 A HIS S 


A 
+ 
Х| 
à 
A 
y 
К 
N 
^ 

















39 427 4 
42 « E 2273 ж; 27 
BS 62-07%, 

de a ae (66) 








[n Late ln jr „жЕ 


he Mes Tyg 





Equations (64), (65), and (66) are the three dimensionless 


equations of motion of the system and contain the following eight 


dimensionless parameters 


x 





№ ^h x 











М 

MM 229 77, ж /73 A 
k: Ar {> 20 
Es 229 2,97 Ж 
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ПІ. NUMERICAL INVESTIGATION OF STABILITY 


Rather than attempt to find a solution for Zquations (64), 


(65), and (66), it was decided to introduce selected values of the 
46 


first six dimensionless parameters and, treating 2 and л 


4 





as variables, to use Routh's Stability Criteria to determine 
the range of combinations of these two variables in which the 
system would be stable. 

To simplify the writing of the determinantal equation 
necessary for the use of the Routh criteria, introduce the fol- 


lowing notation where 


A „ coefficient of Q. in Equation j 


7 
Ж coefficient of Q, in Equation j 
iy coefficient of Q, in Equation j 
and 
j = 1 corresponds to Equation (64) 


j = 2 corresponds to Equation (65) 
j = 3 corresponds to Equation (66) 
With this notation Equations (64), (65), and (66) may now 


be written as 
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с, +в, 0, 48, 0.34: 208 8 751111 
в, + Cu Ë, + B,, Q,+ A, Q,+C,, CO, + By Gt An Ge ee 


B, 0,4 C, G,+B, ©, + A,, Q, + Cys G, + В,, О, +А,, О, = 0 


where, letting 2 = Ж and J = < , the coefficients in Zquation 
2 














(67) аге 

A, = 4/22) 6 “4 

E c ^ 252822 а 

E 5 ES (эж Зо 25 A 

A,, а — (2+ E — — a 

24 = | _/ Hf — M E DEA) = E ta] 
В, = 2 20/ 

Biz = 5 Lale (rs 4) 


B,; ? х A ~ р 


Cys 
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М А Ir E 7 T L^ 252] 





^L di жү" ы. 2% АТ, 232 Ls y 


MI, + Т. 


AI io ze) Ae 
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жж 
(+ 





És 









AS 2,9 


etir (de 
E Al TS 3-1) 


4,7 Er P 


Ф. É, Ma ^fz 
279 АИ a ua с; 
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Inasmuch as oscillatory solutions are of main interest, 


+ 


assume 42 = e? 
Then Фу ге Mei re 


апа a ec DR q. 


For the system of homogeneous differential equations (67) 
to have other than a trivial solution, the determinant of coefficients 


of Q. must vanish. That is 


C, A + В.А B, A+ Ay B,A +A, 
B.A Gud + BENS д Сл А+ В.А + А |=0 (67a) 


Bà С, А+В, х +А,; G.A Eos A, 


expanding the determinant of (67a) 

© ее Y 
Cc, FA +(C, H + B, F)A + (8,H+C,J - B¿S + Bj, W) A 
H(C,K+B, J-B,¿T+B,X)A + (C,R + B, K- B,U+ BY) A 
+ (B,2R - B¿V+B,2)A =0 


which can be written 


РА +рА'+рһ^ +р,х +р А +p- 50 (68) 
where 

р =С„Ё 

р =С„Н+В„Е 

p =B,H+C,J-B,S+B,, W 


p, =C,K+B,J-B,T+B,X 








and 


2 
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B,K- B,U+B,Y 


B.V+B,Z 


Sa Cys 5 CC 


+ Е 


У? 732 


42732 


4432 A, B, 


А,, F A, А, 


Cys + В, В 


ж 33 


y3 


Bayt A Ж; 
А,, Ags P A, Ay, 
Bc - c; 
Х,С, 582 
A,, B, + А, В, 


Ay, Ay = Ay¿ A 542 


ТЭГ? 


«t£ "235 33 
+ В, Б,,- 826, = Arya ч В,, YE 


- A B,,- A, B 


42 9; 797 


- А С,, - В, В 


и 33 393 Уу 


= A, B,, - А,В 


79 W 


-A,C,,- B,B 


34.77 Жу 


-A B,,- A,B 


ж 32 yt  "/ 


The system will be stable if the coefficients of Equation (68) 


and the following terms of Routh's stability criteria are all of the 


same sign (Reference 3) 


l. we 


A045 2 dz autem (69) 
^-^, 


3. 


E 


28 


2227 (BB -RRB XZ 


олду сету 
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NUMERICAL EXAMPLE FOR "V-2" TYPE ROCKET 


To find a numerical solution to the problem the ''V-2"' rocket 
(weight approximately 28,000 lbs.) was chosen as the main rocket, 
| Ма, and a small five-second booster providing an initial accelera- 
tion of 1-1/2 g was selected as the booster, M2. 


The dimensionless parameters of the system were then es- 


timated to be 








^L ғ -< ^ 25 
= Z = LF x — — = - 2, И 
22 иж 27 ж 72, l 
4 зай 26 
RE DO сш x 
2; 27 VB aO 
= O (at launchin 


The stability of the system was then to be investigated for various 


values of the two variables a= 4 and s= A, 
4 A 


Because of the inherent complexity of the constants of Equa- 


tion (69), a vast amount of time would be required to investigate 


& and 4 . How- 


ever, the writer optimistically set out to evaluate these constants 


every likely combination of the variables 


for various combinations of a ands . Inasmuch as no combination 


was found where even the constants p top _ were of the same sign, 
much less where the test functions of Zquation (69) were of the same 


sign, it was decided to reduce the complexity of the system in order 


to expedite the computation. 


2B - 


Accordingly, Q, and а, were assumed to be zero, that is, 
the system was restricted to vertical and rotational movement of 

the center of gravity. | This assumption reduced the system of equa- 
tions to two, namely, Equations (65) and (66), and reduced the labor 
involved in applying the Routh stability criteria substantially. Once 


a stable range of a and 4 had been determined, these values were 


then to have been used in the general case for closer investigation. 
STABILITY OF SYSTEM WITH TWO DEGREES OF FREEDOM 


If both O, and CO, are assumed to be zero, Equations (65) 
and (66) may be re-written (using the abbreviated notation intro- 


duced earlier) as 


са +8,09, + A,,Q, +C d, * B E, +A с, =0 
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C,,@, + B,O, +4,a, +€,0, +B,Q, +4,0, =0 


УР Y 


and assuming Q.«< e" ‚as in the more general case, the determi- 


nant for this case becomes 


£ 
СА + Вл + A, e P E 


2 2 
C, À+ B, A+ Aps C, A B At Ay, 


The expanded determinant may then be written 
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2 3 
-C y A + (Cs. В, + В,С = CUBE. ж B,,Cs2) A 


yz C, 


(быб 


72 43 


ж 
“(С,,А,,%А,С,,%В,В,,- С,,А,,- А,С,,- В.В,)А 


32" «2 


* (8 72 A. + А,В, Е; В.А, = А,В JA * (А,,А уэ” A,¿As7) = 0 (70) 


Furthermore, since the parameter zy equals zero at launching, 
z 


the last term of Equation (70) becomes 


Asa Mas ~%y - =e A AS 22229) 55 2] 


Rs CA у 
Жы (ӨЛУ аа ет. 22 e z aj = 
Therefore Zquation (70) reduces to a third order equation of the 


form 


2462 IAEA E A а (71) 


where 
p *C,C4,-C4,C, 


= C,,B), + B),C,, - Ct, - BEG. 


В. т 22-97 
р <С,4,,%А, Соз % В, В), -С,А,,- А,С,,- В, В,, 
ЛЕӘЗ 5 2285. ЕЕ 


and the Routh stability criteria require that p ,p , P, B,» and 
(8 -4 4) all be of the same sign. (Reference 3). 

The свеїцїёлы of Equation (71) may now be evaluated by 
substituting the parameters previously selected for the "V-2" rocket 


with five second booster. Collecting terms these coefficients may 


be expressed as follows 


fo = Ке) (72) 
where 
О = are 
G = 22802 EL ЛЫ 
P 3 
ET 2 
27- 2/4» <) (73) 
where 
Z = -LPY 077 
Cae — — — г гы 
/ж4 ж 
ГА = GC bad + CZ (74) 
where 
2 - 476 АВ 
4, = 353:10 TD 447310 Y — 444 pa 
@ = - (276 ор Итро za” 73.” 
Ө = еі, бу + Сз (75) 
where 
Ф,- 4024 0 Pa? 
d, - -(& oco th 812 10 53) 
y 
G = 303.103 8.12 PP P e zen 








зэ? 

Before computing numerical values of the coefficients of 
Equation (71) it will be profitable to draw certain conclusions from 
an inspection of Equations (72) through (75). 

a is defined as the quantity ES ‚and g is defined as the 
quantity 4 , where £ andk, are inherently positive quan- 
titles by the physical nature of the koin. Therefore, it can 
be concluded that for any real solution of the physical problem, 
both a and 4 must be of the same sign, that is, either both posi- 
tive, or both negative. 

In Equation (72) the sign of p, is independent of 94 and 
takes the sign of a. Therefore numerical values of p, need not 
be determined at this stage of the investigation since the sign of 
p, will be evident by inspection. у 

In Equation (73) a, is always negative, andc, has roots 
at Æ = 0.456 апі 4 = -1.93, 
Since both à, andc, are 
negative for 4 « 1.93 and 
since a is negative when 7 


is negative, it follows 


that p, is always posi- 





tive in this region. 
Negative a yields negative p, from Equation (72). Therefore p, 


and p, will be of opposite sign for all values of Z < 1.93, and, 
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by Routh's stability criteria, stability is impossible for this case. 
By analogous reasoning stability is impossible in the range 
0<,7< 0.456. 

In Equation (74) a, = 0 when 4,2 3 0, anda, is negative for 
all other values of 4% „ Furthermore, b, has a rootat 4 22.9910 ^, 
апа с, has a root at 
(3 =-2.68x10 . Since 
a, , b, , and c, are all 


negative in the range 


y 





O LA L IY Jo 


it can be concluded that stability is impossible in this range by 
reasoning analogous to that inthe case of p . 

In Equation (75) a, = 0 when 3 = 0, and a, is positive for 
all other values of 8. 
Furthermore b, has a 
root at Z = -1.34x10 
andat Z =0. Also 
C, has roots at /Ф- 
0.456, 3 = -1.92, 


P 





No conclusions can be 
drawn in this case. 


To summarize the conclusions of the preceding paragraphs, 
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the possibility of stability need only be investigated in the ranges 
-1.93</4< 0 and 2.54xl0 4,5 « 0,456. 

In evaluating the coefficients p, , p, , and p, the follow- 
ing procedure was adopted. Selected values of ,Z within the 
positive range 2.54x10 <4 < 0.456 were used to find the con- 
stant coefficients of Zquations (73), (74), and (75). Then the 
positive roots of the resulting equations were computed as func- 
tions of a. Negative roots were not investigated since a must 
be of the same sign as 4% . These values are tabulated in Table 
І. 

Similarly, values of the coefficients and corresponding 
negative roots of Zquations (73), (74), and (75) were computed 
and are tabulated in Table II. 

The data of Tables I and Table II are plotted qualitatively 
in Figure 6 and Figure 7. 

In Figure 8 and Figure 9 are plotted the logarithms of 
maximum and minimum limits of a for which BB. and P, 
are each positive. 

It is apparent from Figure 8 that there is no single value 
of a for which p, , р, , апа р, can all be positive since a,m, 
for positive p, greatly exceeds a,., for positive p, for all 
values of 2.54x10 <,/8< 0.456. 


Likewise, in Figure 9 there is no single value of a for 
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which p , p, , and p can all be negative since 


la,,,| for nega- 
tive p, 


greatly exceeds "ja ,... | for negative p, for all values of 


Therefore, it is clear that quations (65) and (66) do not 


have a stable solution for any combination of a and 4 for the 


parameters chosen. 


STABILITY CF SYSTEM WITH REDUCED MASS RATIO 


In view of the negative results obtained in the previous 


case it was decided to investigate the effect of changing the par- 
mM, 


ameter 7° . The following set of dimensionless parameters 


was then selected 








ШІ 2 222 — ҮЭ th Te = rt 

ж: A3 жж” 

24 2048 — = уво“ 

Ra tir Y 
fe = e a r 

= t launchin 
2” ( 8) 
; € 

Note that 2 was decreased from 7.75 to 1,5, > 


227 
228 was in 
4/9 


creased from 2.48xl0 ^to 8.75x10 . The remaining parameters 


E -< 
was increased from 1.96x10 to6.72x10 ‚апа 
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were unchanged. Thus, this second set of parameters describe 
a system essentially the same as that first investigated except 
that the mass ratio -22- has been markedly reduced. 
The coefficients of Equation (71) are again evaluated by 
substituting these new parameters. Collecting terms the coef- 


ficients may be expressed as follows 


P = ¿(4d +<) (76) 
where a, = sor: /o” 
c = 282.00, A FO vom? 
e a? 
Р, = (46° ey) (77) 
where a, = 777+ 70” 
87 -7 25 
с 2 42486470  /4/74/2Ф  5:64./0 
4 a? Ж 
Р, = ее» б, + СЕ (78) 
= е 2 
where a, = - 422 :/0 4 
bi: 13100 82.480 8 6.34. 077 
2 
с, = — (entm {нё е ч ханыг”, A 
P, = au + bz £ + Сз (79) 


-е „А 
where a= 232 :/0 3 


(MES AA 822.073) 


c 
w 
n 


ж 


б 
n 


-/ж- - ^ - 
232.70 A. 22470 Азр 4 1,77 70 
Fr 
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AS previously stated, a and, must be of the same sign. 
Again, the signof p, is independent of 4 and takes the sign of a. 

In Equation (77) a, 

ee 

is always positive. Also 
c, has roots at ,2 - 0.68 Ж 
and at 7 = -2.84. 

In Equation (78) 
а, = 0 when 7=o, and 
Фа is negative for all 
other values of 7 . 
Furthermore b, has roots 


at 7= -5.6x10° andat 


4 = 2.55x10 7. Alsoc, 





has a root at g= -1.12x10° . 
Since a, , b, , and c, are all negative in the range 0 «,4 < 2.55х10 
stability is not possible in this range. 

In Equation (79) a, = 0 when g = 0, and is positive for 
all other values of, s . Furth- 
ermore b, has roots at,22 0 
and at,4 - -5.62x10” . Also, 


c, has roots at ,Z 3 0.625, 





8 = -2,83 апа 5 = -1.12х10° . 
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It may be concluded that the possibility of stability need 
only be investigated outside the range 0«,7« 2.55x10 7, 

As in the previous case selected values of both positive 
and negative 4 were substituted into Equations (77), (78), and 
(79), to compute the constant coefficienta of Equation (71), and 
the resulting equations were then solved for maximum and min- 
imum allowable values of a. The results are tabulated in Table 
III and Table IV, and are plotted qualitatively in Figure 10 and 
Figure ll. 

Figure 12 and Figure 13 are again logarithmic plots of 
the maximum and minimum limits of a versus the logarithm 
оғ”. 

In Figure 12 it can be seen that stability is possible in 


the range 10:22 « 4.5х107 only if a 


har 


for positive p, is 
greater than a,,,, for positive p, . (The curves appear to co- 
incide.) However, an analysis of Table Ili reveals that for all 
values of 10 7,4 « 0.68 both c, and c, are negative whereas 
a, is negative while a, is positive in this range. Therefore 


Ёх Їнэ ie —F— while a, „р ШЕ 2 | 


Thus it follows that a is greater than a in this range 


I win 2 max 


and stability is impossible. 
Similarly, in Figure 13 stability is only possible if 
ja 


E for negative p, is greater than Ja.) for negative p, . 
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However, an analysis of Table IV reveals that for all values of 


-2.8 « 3< -107 both C, and c, are positive while a, is negative 


and a, is positive in this range. Therefore,|~, , || — /ж БЕЗІ 


while nal 21 Ё үс | and thus |4;,,,| is greater 


than Mames] > 

Therefore it is only necessary to investigate more close- 
ly the range 0.68 <4< 4.5xl0^ and -10 <7 < -2.8. 

Returning to Equation (71) and multiplying the analytic 


expressions for the coefficients of p and p, the following equa- 


tions are derived if only first order terms are retained. 


A = RK tbk ta 
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Comparing Equation (80) with Equation (82) it can be seen 
that botha, andb, are multiples of à, and b, respectively by the 
factor 2 — . For positive, 3 , stability exists only if o, „ma: 


ig greater thana,,, » that is, if 2 < = абас. 
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But the factor 4% “%Zw is numer- 
Ж”. LT, + ^T 


ically negative so that the condition for stability can be written 


Сз 

ES == iion NR 

> Mi Sale 
^7 17,417 





Substituting the values of c, and c, from Equation (80) and Equa- 


tion (82) this condition becomes 


LEO. Mi 
777, 27, Zig 
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which, upon collecting terms, becomes 


27222228 
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Inserting the numerical values of the parameters, the second term 





= -6 
of therdenominator is of the order 10 since Ж = СО 
Therefore this term can be neglected in comparison with the first 


term of the denominator and the expression simplifies to 








> E) (83) 





= 


Solving this inequality with the parameters chosen, the result is 


Ят» 16 








=l- 
^s an example, choose 4 = 2. Inserting values from 
Table Ill into Zquation (78) and Equation (79), for Y = 2 the 


following values of a are obtained 


е ЛИЕ ВО 
— 


a 


се 





до 26.00” (7 =L SY P 


dy x test 2 SS ж 

Thus p, changes from positive to negative when 
a, = 5.6х10”(1-2.54ж19 2, while p, changes from negative to 
positive when a, = 5.6x10"(1-9.85x10%, and both p, and p, are 
positive in the range (5 6x10 -5 .52)<а< (5 6x10" -1, 42) . Re- 
sults with a similar sensitivity are obtained for larger values 
of g but smaller values of a. (It can be seen from Figure 12 
that the product cf « and g is approximately constant and equal 
to 10°.) 

Chere is yet to be applied the final Routh criterion, 
namely that 4 - Е be positive also. Choosing p, at its 
smallest allowable value, say 0,» and calling the correspond- 
ing value of p, sorne small positive value, say £, A- — E, 


Chus all the stability criteria are satisfied and the system is 


stable, Therefore, it can be concluded that the system is stable 








that р, and p, both be negative. Therefore stability is not 


possible for negative values of 4. 
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IV. DISCUSSION CF RESULTS 


Only two possible sets of parameters for a launching sys- 
tem of this type have been investigated. In the first case where 
the mass ratio was 7.75 it was found that stability was impossible, 
whereas, in the second case where the mass ratio was 1.5 it was 
found that stability was theoretically possible for values of 4 be- 
tween 1.62 and 4.5xl0” if the effect of lateral motion of the center 
of gravity of the system was neglected. For any given value of 4 
in this range stability was possible for only one correspnnding 
value of a. 

Returning to the physical aspects of the problem consider 
now the meaning of the variables a and s . 4 is defined as E 
where 4 is the distance from the center of gravity of the booster, 
M,, to the pin connecting it to the strut to the main rocket, M, , 
and k, is the radius of gyration of the booster, M, . For any sys- 
tem of this type k, will сан be of the order of one foot or 
greater. Therefore, 4 =4 4, will always be equal to or greater 
than s. That is, as a first approximation, 4 is approximately 
numerically equal to Z. Thus, from a practical consideration 

4, , and therefore Z , must certainly not exceed some small num- 
ber of the order of 1 to 10, depending on the size of the booster, 


for otherwise the pin will be located too near the booster exhaust, 
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Now a is defined as $ where £ is the length of the strut 
from the pin in the booster M, to the center of gravity of the 
main rocket, M,. Here again there is a practical upper limit 
for £ , first, because the strut was assumed to be of zero mass 
in the derivation of the equations of motion, and second, because 
the strut was assumed to be of infinite rigidity. A long strut 
would present structural problems in rigidity since it is essen- 
tially a column acting under both axial and bending loads. Re- 
stricting A. to values of the order of a few feet it is clear that a 
should not exceed 100 or 200 at most, since these values corres- 
pond to an -% of 25 to 100 feet. 

In view of the foregoing remarks consider the product of 
the two variables a and 4. Thus 27 -44 = E ‚ ик, 
is of the order of 1 to 10 feet, say, and 4 is restricted to 25 to 
100 feet, then it can be seen that the product a# should not ex- 
ceed 25. Returning to Figure 12 note that, in the range 1.62<4< 
4.54x10” , a4 is approximately constant and equal to about 10 2 


Thus, even though theoretically stable, this system is of no prac- 


tical significance for the mass ratio investigated. 
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V. CONCLUSIONS 


1. A compound pendulum system of launching stabilization is not 


stable for a mass ratio 7.75. 


2. Such a system is theoretically stable for a mass ratio 1.5 if 
the effect of lateral motion is neglected. However, even in this 
restricted case the system is of no practical significance due to 
the excessive length requirements for the strut from the main 


rocket to the booster. 


3. Since the system is quite sensitive to the particular param- 
eters chosen, several more investigations of specific cases are 
required before any general conclusions can be drawn as to the 


practicality of this method of stabilization. 
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Figure 6. Variation of p , p, , and p, with a for several values 


of positive 8. E = 7.73) 
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Figure 7. Variation of p , Р, >» and р, with a for several values 


of negative 4. ( E 57.15) 
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Figure 8. Minimum and maximum limits of a which yield posi- 


tive p. for positive g . ( — = 775) 
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Figure 9. Minimum and maximum limits of a which yield nega- 


tive p for negative 4. ( == 7.75) 
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Figure 10. Variation of p , p, , and p, with a for several values 


of positive 4 . (zs 1255 
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Figure ll, Variation of p , p ,andp with a for several 


values of negative 8. ( ==1.5) 























Figure 12. 
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Minimum and maximum limits of a which yield positive 


р, for positive 4. ( ж 21.5) 
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Figure 13, Minimum and maximum limits of a which yield nega- 


tive p for negative 4. (4 =1.5) 
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Figure 14. Schematic diagram of system showing inertia 
forces for use in Newtonian derivation of equa- 


tions of motion. 
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APPENDIX 


NEWTONIAN DERIVATION CF EQUATIONS CF MOTION 


The equations of motion for the two dimensional system 
may be derived by applying Newton's Second Law of Motion pro- 
vided that proper consideration is given to the inertia forces of 
the system. In Figure 14 these forces are indicated by broken 
lines. In this derivation the main rocket, M, , and the booster, 
M; , will be treated separately as free bodies and the reaction 
of the pin connection will be replaced in Figure 14 by its two 
components, f, in the vertical direction, and f, in the horizon- 
tal direction. The pin is assumed to be frictionless and there- 
fore no moment is transmitted. 

Treating the main rocket, M,, as a free body first, the 


sum of the vertical forces acting on it are 
ZV= / -Mg-Mg 
and therefore 


«7 7 | 0) 


Similarly, the sum of the horizontal forces acting on M, are 
ZH nS 
and therefore 


4 = аў, 
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Finally, the sum of the moments acting on M, are 


EM= ZA +h Moose -4Cs0g, 


Ж- LA sin ga - 4.6 ғз 2 (3) 


Treating the booster, M,, as a free body, the sum of the 
vertical forces acting on it are 

IV = “соғу, Brin E A AAA Y = 22 
and therefore 

— Feos g, - 2 sin gu -1LF -1h fy = JA (4) 
Furthermore, the sum of the horizontal forces acting on M, are 

ТН. Fang + Dees fy Aa 27222 
and therefore 

Ж Fain M ELT (5) 
Finally, the sum of the moments acting on M, are 

гм=& 6 + „ў, РА 207 AA cos 7 7 2/5 5) 
or 

4% "LZ - 2-6 ез -4-6 «ту; - 2/4-.6) (6) 

Substituting the results of Zquation (1) into Equation (4) to 
eliminate f, 

Mg tif, ео ў Рел Ig MR I, (7) 
Similarly, substituting the results of Equation (2) into Equation (5) 


to eliminate f, 


Mg, < эт) Рез u a (8) 





Bi- 
Substituting the results of Zquations (4) and (5) into Equation (3) 


to eliminate f, and f, 


1% = [Fcos Р зоол 149 fa fe) э” Ж 


(9) 
(mg + Leos — > -1% fp) £ 44525 
Similarly, substituting the results of Equations (4) and (5) into 
Equation (6) to eliminate f, and f, 
4% * 4% (Fur Des ta af) er 
[Feos Ф -Psa fo - 1% 7-14 Ja) mn fy 
3 204 tG) (10) 


Assuming the angles q, and q; sufficiently small such that 
terms of order q* may be neglected, the following approximations 
can be made | 

23% > 7. апа 992 = / 
Therefore, Equations (7), (8), (9), and (10) may be rewritten in 
simpler form. Thus 
From Equation (7) 

ж EE - (7 A (11) 
From Zquation (8) 

КАА с. aa 


From Equation (9) 


273 = [F -Lpy UGR ARE 


ee (13) 
Ta “2-2 2 - 2 fa) A 
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From Equation (10) 
LZ + ER = (ar A A 
—-(“- 29 -EI -“% Aja, -A E) (14) 


The following approximate relations between the principal 
coordinates q , а, , 9, , and q, , and the extraneous coordinates 


а, , Ч, 4 а, 4 Ч. were proved earlier in Part I of this paper 











en » 
т СС (16) 
С ЕЯ - 3 [(4%,- A) (17) 
jo * fi - a, (4-4) (18) 


Taking the first and second time derivatives of Equation (15) 


Ў = Л д8)” AR 45) (19) 


and 
i = gap) pCR E 
— — 44) % хо lt, (“Gs — 4) (20) 
Taking the first and second time derivatives of Zquation (16) 
Ё бе аа = 
and 


Ё 2 Б зурган 1-4) (22) 


Taking the first and second time derivatives of Zquation (17) 
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$ = Е "eA A fu) - ж mom (^ 5 7 AR) (23) 
and 
BR At IA oiii AL 3) 
“аг Aja) - OZ - 47, ) (24) 
Taking the first and second time derivatives of Jquation (18) 
A (25) 
and 


> - 2 ж/е Es 

A (26) 

Equation (11) can now be rewritten in terms of the principal 
coordinates only by substituting the results of Equations (22), (25) 


and (26) for the extraneous coordinates 


"е E ааа /4- - Яр - <-90%%4)-2р, (27) 
Note that Zquation (27) is identical with Equation (27) of Part II of 
this paper. 

In a similar manner, Equation (12) can be rewritten in terms 
of the principal coordinates by substituting the results of Equations 


(20), (23), and (24) for the extraneous coordinates 


(28) 





>” 7, 3 е >< 20: 
au A P 


which is identical with Equation (23) of Part II of this paper. 
Substitution of Equations (23), (24), (25), and (26) in Equa- 


tion (13) yields 
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- 4 jn) * яо СОЗ” 4 fe) - ge — x f 








G2 


4, (4-4) 





МДА — 
01-7) ар-ар) Go fe + 22 AA 


= "zz LA + 2/ 


^A, 


(29) 


which is identical with Equation (44) of Part II of this paper. 


Finally, substitution of Equations (23), (24), (25), and (26) 


in. Equation (14) yields 





Д ADA) A A OB A) 


ел y Lr, M Ugo) A RÀ 


Y OM 
20-24) х - D l + Ns 








(30) 





is Ap 


which is identical with Zquation (59) of Part II of this paper. 
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- Acceleration of gravity 

- Moment of inertia 

- Specific Impulse 

- Radius of gyration 

- Mass 

- Mass rate of flow 

- Lagrangian generalized space coordinate 

- Dimensionless Lagrangian generalized space coordinate 
- Lagrangian generalized force 


- Kinetic energy 
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